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Abstract 

With the first measurements of the branching ratios and the direct CP asymmetries 
of Bg — > K~K^, K^TT^ decays by the CDF collaboration, we constrain the relevant 
parameter space of the Minimal Supersymmetric Standard Model with R-parity violation. 
Using the constrained R-parity violating couplings, we further examine their possible 
effects in Bg K~*tt^, K^*^~ p'^ and Er(*)='=/f {*)t decays. We find that some branching 
ratios and CP asymmetries are very sensitive to the R-parity violating couplings. The 
direct longitudinal CP asymmetries of tree-dominated process Bg K*~ p'^ could be 
enlarged to ^ 70% and the longitudinal polarizations of Bg — > K*~K*^, K*~ p^ decays 
could be suppressed very much by the squark exchange couplings. Near future experiments 
at CERN LHC can test these predictions and shrink/reveal the parameter spaces of RPV 
SUSY. 



PACS Numbers: 12.60.Jv, 12.15.Ji, 12.38.Bx, 13.25. Hw 

* Corresponding author, yangyd@iopp.ccnu.edu.cn 
t Mailing address 



1 



1 Introduction 



In the recent ten years, the successful running of B factories BABAR and Belle has provided rich 
experimental data for and B^, which has confirmed the Kobayashi-Maskawa CP asymmetry 
mechanism in the Standard Model (SM) and also shown hints for new physics (NP). Among the 
rich phenomena of B decays, the two-body charmless decays are the known effective probes of 
the CP violation in the SM and are sensitive to potential NP scenarios beyond the SM. The two 
body charmless Bg decays will play the same important role for studying the CP asymmetries 
(CPA), determining CKM matrix elements and constraining/seraching for the indirect effects 
of various NP scenarios. 

Recently the CDF collaboration at Fermilab Tevatron has made the first measurement of 
charmless two-body Bs decays 

B{Bs K'K^) = (24.4 ± 1.4 ± 3.5) x 10"^ 
B{Bs A'~7r+) = (5.0 ± 0.7 ± 0.8) x 10~^ 

Acp{Bs K'Ti+) = 0.39 ± 0.15 ± 0.08. (1) 

The measurement is an important mark of Bg physics, and also implies that many Bg decay 
modes could be precisely measured at the coming LHCb. 

Compared with the theoretical predictions for these quantities in Refs. 

a. a 

and 



based on the QCD factorization approach (QCDF) ^, the perturbative QCD (PQCD) jlO|. 
and the soft-coUinear effective theory (SCET) HI], respectively, one would find the experimental 
measurements of branching ratios agree with the SM predictions within their large theoretical 
uncertainties. However, NP effects would be still possible to render other observable deviated 
from the SM expectation with the branching ratios nearly unaltered |l^ . 

The related decays B^ — >■ K^*^~ K^*^^ , K^*^~7i^, K^*^~ have also been extensively studied 



in the literature 



13 



14 



151 ]. The four decays Bg are governed by 

the b suu transition at the quark level, which are penguin-dominated processes. The tree- 
dominated decays Bg K^*'^~7r~^, K^*'^~p~^ are induced by 6 ^ nn J where the direct CPA are 
expected to be small in the SM. At present, among many measurements of B^ a decays, several 
discrepancies with the SM predictions have appeared in the corresponding penguin-dominated 
b — > sqq {q = u, d, s) processes and tree-dominated b — > dq'q' {q' = m, d) processes. For example, 



B ^ TiTT, 71 K puzzles 
B pK*,(f)K* decays 



13, 
22 



19 



201 ] and the large transverse polarization anomaly in 



231]. Although the discrepancies are not statistica' 



there is an unifying similarity pointing to NP (for example, 12 



24 



25 



26 



27 



ly significant, 



29, 



30|). 



There could be also potential NP contributions in B^ K^*^ K^*^~^, K^*^ 7r+, K^*^ p+ decays, 



which have been analyzed with different NP models p^, 130|, 131|, |32| . The measurements given 
in Eq. ([T]) will afford an opportunity to constrain NP scenarios beyond the SM. 

Among the NP models that survived electroweak data, one of the respectable options is the 
R-parity violating (RPV) supersymmetry (SUSY). The possible appearance of the RPV cou- 
plings 33|, l3J] , which will violate the lepton and baryon number conservation, has gained full at- 



tentions in searching for SUSY [35|, l36|, l37|, l38| . In this work, we will study the Bg ^ 
/^(*)-7i-+ and K^*^~ decays in the Minimal Supersymmetric Standard Model (MSSM) with 
R-parity violation by employing the QCDF. The four Bs x^*^^ K^*^~ decays are all induced 
at the quark level by 6 ^ suu process, they involve the same set of RPV coupling constants. 
The Bg A'(*^~7r+, K^*^~ decays are due to 6 duu at the quark level, and they also involve 
the same set of RPV coupling constants. Using the latest experimental data and the theoretical 
parameters, we have derived new bounds on the relevant R-parity violating couplings, which 
are consistent with the bounds from Bu,d decays. With the constrained parameter spaces, we 
predict the RPV effects on the other quantities in B^ ^i*)- ^ K^*'>~'k^ and A'(*)~p+ 
decays which have not been measured yet. We find that the R-parity violating effects on some 
branching ratios and direct CPA could be large. For example, the squark exchange couplings 
could enhance the direct CP asymmetry in the longitudinal polarized mode of B^ K*~ p^ to 
~ 73% and suppress the longitudinal polarization fractions of Bg K*~K*^ and K*^ p^ to 
~ 0.5. The mixing-induced CPA are also found to be sensitive to the RPV effects. Therefore, 
with the ongoing B-physics at Tevatron, in particular with the onset of the LHC-b experiment, 
we expect a wealth of Bg decays data and measurements of these observables could restrict or 
reveal the NP parameter spaces in the near future. 

The paper is arranged as follows. In Sec. 2, we give the expression of the CP averaged 
branching ratios, the direct CPA, the mixing-induced CPA and the polarization fractions within 
the QCDF approach in Bs A'(*)-A'(*)+, A'(*)-7r+, A:(*)-p+ systems, where the RPV SUSY 
effects are included. We also tabulate the theoretical inputs in this section. Sec. 3 deals with 
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the numerical results. We display the constrained parameter spaces which satisfy the present 
experimental data of Bg decays, and then we use the constrained parameter spaces to predict 
the RPV effects on the other observable quantities, which have not been measured yet in 
B, K'-*^-7i+ and K'^*^- p+ decays. Sec. 4 contains our summary and conclusion. 



2 The theoretical frame for B M1M2 decays 
2.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the AB = 1 transition at the scale n ~ rrib 



is given by 



39| 



^ / 10 \ 

'^eff = ^ E C'lQ? + C2QI + E C^Q^ + Cj.Qj, + CsgQsg +h.C., (2) 

p=u,c \ j=3 / 

here = VpbV*g hi b q transition {p G {m, c}, g G {d,s}) and the detailed definition of the 
operator base can be found in 3^. 

With the weak effective Hamiltonian given by Eq. ([2]), one can write the decay amplitudes 
for the general two-body hadronic B — M1M2 decays as 

A^^{B ^ M1M2) = (MiM2|Hf//|B) 

= EE^pC'.(/i)(MiM2|g.(/i)|i?). (3) 

P i 

The essential theoretical difficulty for obtaining the decay amplitude arises from the evalua- 
tion of hadronic matrix elements {MiM2\Qi{^)\B) , for which we will employ the QCDF 
throughout thispaper. 

The QCDF [9| allows us to compute the non-factorizable corrections to the hadronic matrix 
elements {MiM2\Qi\B) in the heavy quark limit. The factorization formula reads 

(MiM2|g,|5) = {Ff-''% * fM,<^M, + [Ml ^ M2]) + Tl' * fB<^B * M<f A/1 * /m.$m„ (4) 

where p^^'^^^ is the appropriate form factor, are leading-twist light-cone distribution am- 
plitudes and the star products imply an integration over light-cone momentum fractions. By 
the above factorization formula, the complicated hadronic matrix elements of four-quark oper- 
ators are reduced to simpler non-perturbative quantities and calculable hard-scattering kernels 

rpI,II 



4 



Then the decay amphtude has the form 

j^SM^Q ^ M1M2) = A4anM2|J2|0)(Mi|Ji|B) + 6?(MiM2|J2|0)(0|Ji|fi)|, (5) 



p 

P i 



where the effective parameters af including nonfactorizable corrections at order of a^. They 
are calculated from the vertex corrections, the hard spectator scattering, and the QCD penguin 
contributions. The parameters 6f are calculated from the weak annihilation contributions. The 
factorized matrix element is given by 

Am,m, = (M2|(g27M(l-75)g3)|0)(Mi|(67'^(l-75)gi)|S), (6) 

which can be expressed in terms of the corresponding decay constants and form factors. We 
will use the QCDF amplitudes of these decays derived in the comprehensive papers g], Q as 
inputs for the SM amplitudes. 

2.2 R-parity violating SUSY effects in the decays 



In the most general superpotential of MSSM, the RPV superpotential is given by 



33| 



yVnpv = ^iUHu + -\[ij]kLiLjEl + X'ij^.LiQjDl + -A-y^jf/fD^Z)^, (7) 

where L and Q are the SU(2)-doublet lepton and quark superfields and E'^, U'^ and D"^ are the 
singlet superfields, while i, j and k are generation indices and c denotes a charge conjugate 
field. 

The bilinear RPV superpotential terms fiiLiHu can be rotated away by suitable redefining 



the lepton and Higgs superfields (35j. However, the rotation will generate a soft SUSY breaking 
bilinear term which would affect our calculation through penguin level. However, the processes 
discussed in this paper could be induced by tree-level RPV couplings, so that we would neglect 
sub-leading RPV penguin contributions in this study. 

The A and A' couplings in Eq. (171) break the lepton number, while the A" couplings break the 
baryon number. There are 27 A'j^-^ couplings, 9 Xijk and 9 A^'^-^, couplings. \[ij]k are antisymmetric 
with respect to their first two indices, and ^-yk] antisymmetric with j and k. 

From Eq. ([7]), we can obtain the relevant RPV effective Hamiltonian as shown in Fig. [T] 






Figure 1: RPV tree level contributions to 6 — > uuq process. 
A" A"* 

- [{dk-i^PRUj)i{ua^,PRdi)i - {dk-i^PRU,j)^{ua^PRdi)^ } , (8) 

where Eq. ([H]), Pl = ^~2^^ Pr = V = 'a^imi] /^o = H — f'^/- The subscript of 

the currents (j^)i,8 represents the current in the color singlet and octet, respectively. The 
coefficients ?^~^/^o and r^"^/^" are due to the running from the sfermion mass scale (100 GeV 
assumed) down to the scale. Since it is always assumed in phenomenology numerical display 
that only one sfermion contributes at one time, we neglect the mixing between the operators 
when we use the renormalization group equation to run Ti-^fj^ down to the low scale. 
The RPV amplitude for the decays can be written as 

^R^^(5^MiM2) = (MiM2\nff/\B). (9) 

Generally, the product RPV couplings can be complex and their phases may induce new con- 
tribution to CP violation, which we write as 

AijkALn = \A^JkAlmn\ c"^^^^ , Kjk^lmn = \A^JkAlmn\ e"*^«^^. (10) 

The RPV coupling constant A G {A, A', A"}, and (pRpv is the RPV weak phase, which may take 
any value between — vr and vr. 

For simplicity we only consider the vertex corrections and the hard spectator scattering in 
the RPV decay amplitudes. We ignore the RPV penguin contributions, which are expected 
to be small even compared to the SM penguin amplitudes, this follows from the smallness of 



the relevant RPV couplings compared to the SM gauge couplings. __T 
RPV couplings are insensitive to the inclusion of the RPV penguins 



lus, the bounds on the 



40[ . We also neglected the 



annihilation contributions in the RPV amplitudes. After Fierz transformations, the relevant NP 
operators due to squark exchanges are (^7^(1 +75)^) (67^(1 +75)^) and (^7^(1 — 75)^) (67^(1 + 

6 



75)5')- The factorized matrix element of these new RPV operators is given as follows 

^Wm. = (M2|(n7^(l+75)g)|0)(Mi|(67ni + 75)^)|5), (11) 



-I < 







if Ml 


= M2 = P, 








if Ml 


= V, M2 = P, 








if Ml 


= P, M2 = V, 








if Ml 


= M2 = V and h = 


0, 


mBrriMzfAh 




if Ml 


= M2 = V and h = 


±, 



(12) 



with 

Pf-^Hg^) = fl + ^) Af-^'^iq') tU-^) V^-'^^iq'). (13) 
V ms / V 1713/ 

Using QCDF, we can obtain the RPV amplitudes of P, K^*)- K^*)+ ^ i^(*)-7r+, K^*^- p+ 

decays. There are two independent RPV amplitudes, given by 

\//* \// \'* \' 

A^'^{B. K-K-) = -^f^r^/^i^.x^W - ^V-'^'^r-^N'iK)A',^, (14) 

d ^Li 

^^^^(P.-Ji+vr-) = -^^^f^V"'^''F^^A'^.^-^^^ (15) 

ami 8m| ^ 

with iV'(M2) = 1 if M2 = P and iV'(M2) = if M2 = V". PmiMj is defined as 



V'(M2) + ^H\M,M2) , (16) 



where l^'(M2) and H'{MiM2) are the one-loop vertex corrections and hard spectator inter- 
actions for the new RPV operators, respectively. The RPV amplitudes for P^ K~^K*~, 
Bs K*+K- and P, K*+K*- are obtained from Eq. by replacing (KK) ^ (KK*), 
(KK) {K*K) and (KK) {K*K*), respectively. The RPV amplitudes for P, K+ p*- , 
Ps ^ K*+n~ and P^ ^ K*^p*~ are obtained from Eq. ([IS]) by replacing {Kn) (Kp*), 
(Kir) {K*7i) and (Kn) {K*p), respectively. 

As for V^'(M2) and ^'(MiMs) for M1M2 = PP,PV,VP,VV cases, the exphcit results are 
same as these of SM operator ("^7^(1 — 75)5) (67^^(1 —75)^) except ones for P VV and h = ± 
case. And we get 

V'-{V) = 0, V'+{V) = V^{V). (17) 
H'-{VV) = 0, H'+{VV) = -H'{VV). (18) 
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2.3 The total decay amplitude 

From the above discussions, the total decay amphtude are then given as 

A{B, ^ M1M2) = A'^^'iB, M1M2) + ^^^^(5, ^ M1M2). (19) 
The corresponding branching ratios read 



B{Bs ^ M1M2) = \AiB, ^ M,M2)f , (20) 



where tb^ is the Bg hfetime, \pc\ is the center of mass momentum in the center of mass frame of 
the Bg meson. In the B — > VV decay, the two vector mesons have the same hehcity, therefore 
three different polarization states are possible, one longitudinal and two transverse, and we 
define the corresponding amplitudes as Ao,±. Transverse (^||,±) and helicity {A±) amplitudes 
are related by A\\^± = "^^^^ ■ Then we have 

\A{B^vv)f = iAi'+iAr+i^-r = iAi'+i^iii'+i^±r. (21) 

The longitudinal polarization fraction Jl is defined by 

hiB^VV) = ^ = , , „ , 1^°'' , , , „ ■ (22) 



For the CPA of neutral B meson decays, there is an additional complication due to B^ — B^ 
mixing. There are four cases that one encounters for neutral B^ decays, as discussed in Ref. 



41 



42 



43 



44|. 



Case (i): B^ ^ f, B^ ^ /, where / or / is not a common final state of B^ and B^, for 
example 5° K~ tx+ , p+ , K*^ n+ , K*~ p+ . 

Case (ii): B^ ^ {f = f) ^ B^ with f^^ = ±/, involving final states which are CP 
eigenstates, i.e., decays such as K^K^. 

Case (iii): B^ ^ {f = f) ^ B^ with f^^ ^ ±/, involving final states which are not 
CP eigenstates. They include decays such as B^ — > {VV)^, as the VV states are not CP 
eigenstates. 

Case (iv): B^ if^f) ^ B^ with Z*^"^ ^ /, i.e., both / and / are common final states 
of 5° and 5°, but they are not CP eigenstates. Decays fi°(fi°) ^ K*-K+,K-K*+ 
belong to this case. 



For CP case (i) decays, there is only direct CPA A'^'^p since no mixing is involved for these 
decays. For cases (ii) and (iii), their CPA would involve — mixing. The direct CPA Ac^p 
and the mixing-induced CPA are defined aJ3 

Ak4ir,-nQ r\ \'^k\ ~ 1 ,k,mixfTjO , s\ 2Im(Afc) , , 

Acp [B Jj = ,2 , , > Acp [B ^ f) = [26) 
Afe + i Afe + i 



where k = 0, ||,_L for i? VV decays and k = ioi B ^ PP,PV decays, in addition, 
Xk = l il'^oZfi for CP case (i) and = f ^^{fe^ for CP cases (ii) and (iii). 

;he four time- 



Case (iv) also involves mixing but requires additional formulas. Here one needs t 
dependent decay widths for B^{t) /, B^t) /, B°{t) f and B\t) f 



These time-dependent widths can be expressed by four basic matrix elements 4^ 



42 



43 



4^. 



g = iflUefflB'), h = {f\neff\B'), 

g = h={f\neff\B'), (24) 

which determine the decay matrix elements of B^ fhf and B^ fhf at t = 0. We will 
also study the following quantities 

(S°&5° ^ /) = ^f^, (S°&E° ^ /) = (25) 

A'cf\B'kB' - /) = ^Ipf^, A'c7\B'kB' f) = (26) 

|A^| + i |A^| + i 

with A';, = l{h/g) and A'^' = lig/h). The signature of CP violation is V{B\t) -> /) ^ 
T{B^{t) f) and V{B\t) ^ f) ^ V{B^{t) /), which means that A^Sp'\B^^B^ /) 
-^^f (50&50 ^ /) and/or ^^™"(50&fiO ^ f) ^ -Al}p"{B'>kB' ^ /). 

2.4 Input Parameters 

The input parameters are collected in Table I. In our numerical results, we will use the input 
parameters which are varied randomly within la range. The Wilson coefficients Cj are eval- 



uated at scales fi = rrib (39|]. For hard spectator scattering, we take fih = ^JAQCD'mb■ When 
we study the RPV effects, we consider only one RPV coupling product to contribute at one 
time, neglecting the interferences between different RPV coupling products, but keeping their 



^ We use a similar sign convention to that of [45| for self-tagging i?" and charged B decays 
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interferences with the SM amphtude. We assume that the masses of the sfermions are 100 GeV. 
For other values of the sfermion masses, the bounds on the couphngs derived in this paper can 
be easily obtained by scaling them by factor p = ( ioo"gcv )^- 

3 Numerical results and Analysis 

Now we are ready to present our numerical results and analysis. First, we will show our es- 
timations in the SM with the parameters listed in Table I and compare with the relevant 
experimental data. Then, we will consider the RPV effects and constrain the relevant RPV 
couplings from the experimental data. Using the constrained parameter spaces, we will give 

Table I: Default values of the input parameters and the itlo" error ranges for the sensitive parameters 
used in our numerical calculations. 





= 5.366 GeV, m , = 0.892 GeV, m , = 0.494 GeV, 






m , 


= 0.140GeV, rup = 0.775 GeV, mb(m;,) = (4.20 ± 0.07) GeV, 






m„(2GeV) = (0.0015 ~ 0.003) GeV, mrf(2GeV) = (0.003 ~ 0.007) GeV, 






m,(2GeV) = (0.095 ± 0.025) GeV, r^^ = (1.530 ± 0.009) ps, r^^ = (1.4371 


-0.030N 
0.03lJ P^- 


[46] 


\Vud\ 


= 0.97430 ± 0.00019, |K.| = 0.225211°:™, = 0.003441°;™^ 








= 0.225081™, = 0.973501°:™, = 0.040451°:°^, 






\Vtd\ 


= 0.0084ll°:S, IV^.I = 0.03972l°:K, = 0.999176l°:°K, 






a = 


(90.711:^)°, /3=(21.71J:°)°, 7=(67.6li^)°. 




m 


fK = 


: 0.160 GeV, /k' = (0.217 ±0.005) GeV, = (0.156 ± 0.010) GeV, 






h = 


0.131 GeV, fp = (0.205 ±0.009) GeV, = (0.147 ± 0.010) GeV, 






^0 


^^*(0) = 0.360 ± 0.034, ^^*(0) = 0.233 ± 0.022, Af^^^*(0) = O.l; 


31 ±0.025, 




yBs- 


■^^'(0) = 0.311 ± 0.026, Fo^=^-^(0) = 0.30ll]:[jt 




[48,49] 


fSs = 


= (0.245 ± 0.025) GeV. 




m 


\b = 


: (0.46 ±0.11) GeV. 






al = 


0, a^ = 0.20 ± 0.15, a{ = 0, = 0.1 ± 0.2, 






< = 


= 0.2 ± 0.2, af = 0.1 ± 0.3, af = 0.06 ± 0.06, af = 0.1 ± 0.2. 




[6, 13J 
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the RPV SUSY predictions for the branching ratios, the CP asymmetries and the longitudi- 
nal polarization fractions, which have not been measured yet in Bs — > K^*^^ K^*^~^ , K^*^~ti^, 
K^*^~ p'^ systems. 

For CP case (i) decays Bg K^*^"^^ and K^*^^ p~^, we will study the CP averaged branching 
ratios (B), A^p and the longitudinal polarization fractions (/l)- For CP cases (ii), (iii) and 
(iv) decays Bs K^*)-K^*)+^ we will also study A'gf besides B, Afp and /l. For CPA of 
Bs K*~ K*+ , K*' p+ , we only study the longitudinal direct CPA [Acp'') and longitudinal 
mixing-induced CPA (^^'™^). The numerical results in the SM are presented in Table II. The 
detailed error estimates corresDonding to the different types of theoretical uncertainties have 

nn 

been already studied in Refs. p, ll3||, and our SM results of B, A^p and fi are consistent with 
the ones in Refs. flQ. 

Table II: The SM predictions for B (in units of 10"^), A'^'p, and A'Sf in Bs 
K-K+, K*-K+, K~K*+, K~TT+, K*^Tr+, K~p+ decays within QCDF. ^ K*^K+ 

denotes that 5° and 5^ decay to the same final state K*~ . 



Decay modes 


B 


A dir 


Ami'x 


Bs K-K+ 


[0.89, 4.45] 


[0.02, 0.06] 


[0.21, 0.43] 


Bs K*-K+ 


[0.22, 2.14] 


[-0.07, 0.02] 




Bs K~K*+ 


[0.21, 0.65] 


[0.03, 0.10] 




BskBs K*^K+ 




[-0.72, 0.34] 


[-0.30, 0.04] 


BskBs ^ K-K*+ 




[-0.31, 0.73] 


[-0.34, 0.12] 


Bs K~7i+ 


[0.61, 1.47] 


[-0.09, -0.05] 




Bs K*^n+ 


[0.84, 1.72] 


[0.00, 0.02] 




Bs ^ K-p+ 


[1.33, 3.51] 


[-0.02, -0.01] 





Our results of i? ^ PP and PV are obtained excluding the uncertainties of power correc- 
tions parameterized by the quantities Xa and Xh- In the QCDF, the endpoint divergent 
integrals appear in the hard-scattering contributions and in the weak annihilation contri- 
butions, which are treated with mo del- dependent parameters js] Xh = (1 + f)//e*'^^)ln^ 
and Xa = (1 + ^)yie*'''^)ln^, respectively. The different Xa values are allowed for the 
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Table III: The SM predictions for B (in units of 10 Acp^, A^'p^^ and fi in B. 
K*~K*+, K*~p+ decays within QCDF. 



Decay modes 


B 


A L,dir 


A L,mix 
J^CP 


h 


Bs K*~K*+ 


[0.39, 1.71] 


[-0.04, 0.19] 


[0.70, 0.93] 


[0.38, 0.89] 


Bs ^ K*-p+ 


[1.03, 6.23] 


[-0.06, -0.03] 




[0.86, 0.97] 



four cases PP, PV, VP and VV. Our results of i? ^ PP, PV are obtained without the 
uncertainties of power corrections and set Qa = Qh = 0. For two vector final-state meson 
decays, in order to be consistent with the longitudinal polarization fractions around 0.5 
in the penguin-dominated decays B (pK*^ and p~^K*^, maximal annihilation contribu- 



tion are considered {qa = 0.6 and (pA = —40°) in Ref. [l3|. We also consider the large 



annihilation contribution and suggest qh = 0, Qa = 0.6 ± 0.2 and (fA = (—40 ± 10)°. 
The annihilation topology obviously contributes to B, Aqp^ and Aq^^^ besides fi in 
Bs K*~K*+ decay. For example, AcpiB'^s ^ K*~K*+) receives much larger anni- 
hilation contribution than Aqp{B^ K~ K^) does. It is also noted that annihilation 
contribution could cancel voluminous penguin contribution in (^5 ^ K*-K*+). 

For CP case (iv) Bs K*~K^ decay, the final state can come both from a pure Bg and 
a pure the amplitudes for the direct decay Bs K*~ and the mixing-induced 
sequence Bs ^ Bs ^ K*- K+ . We obtain ^^*^(E,&5, K*- K+) ^ -Afp{BskBs 
K-K*+), however, ^™)?(5,&i?, K*-K+) -A^S^^BskBs K-K*+), which imply 
that its direct CP violation is very small, nevertheless its CP violating effect can appear 
through the interference of the direct decay Bg K*~K^ and the mixing-induced decay 
Bs ^ Bs ^ K*^ . In addition, the theoretical predictions for above CP asymmetry 
parameters suffer large uncertainties, which are dominated by the uncertainties of mass 
and the Gegenbauer moments in the expansion of the light-cone distribution amplitudes, 
and also due to the uncertainties of the form factors and the CKM matrix elements. 

In penguin-dominated decay K*^K*~ , as transverse and longitudinal contributions 

can be be of the similar magnitude, the CP asymmetry and the polarization fractions pre- 
dictions will suffer large uncertainties. For example, compared to Acp{B^ K~K~^), 
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AqP^{B^ —>■ K*~ K*~^) suffers quite large uncertainties, which mostly come from the un- 
certainties of the relevant form factors and the weak annihilation parameter Xa- fiiB^ 
K*^K*~^) has quite large allowed range for the same reason as 

• AcT^'iBl K*~K*+) is much larger than ^ K-K+) in the SM. Large differ- 
ence between them arises from chirally-enhanced terms, which give large contribution to 
penguin-dominated decay modes with pseudoscalar final-states. 

• For the color-allowed tree-dominated decays Bg K^tc'^, K*~tc^, K^p~^, and K*~ , 
power corrections have limited impact, and the main sources of theoretical uncertainties 
in the branching ratio are CKM matrix elements and form factors. Their A'^^p and A^'p"^ 
can be predicted quite precisely, and found to be very small (~ 10~^) due to small penguin 
amplitudes. The uncertainty of K*~ p^) is mostly due to the uncertainties of 
form factors. 

Now we turn to the RPV effects in 5, K^*)-K^*)+^ K^*^-7i+, and K^*^' p+ decays. There 
are two RPV coupling products, X'isiX'ui and A-*3A-i2 contributing to four Bg k'^*)- k'^*)+ 
modes, which involve the quark level process h uus. Four decays Bs K^*^^n'^, K^*^~ p~^ 
are due to 6 ^ uud at the quark level, and the relevant RPV coupling products are A']'32A^\2 
and ya^yni- We use the experimental results shown in Eq. ([T]) to constrain the relevant RPV 
parameters. 

Our bounds on A'/3^A']'2i and X'asKu are demonstrated in Fig. [2] (a-b) by using the ex- 
perimental measurement of B{Bs K"K'^) within Icr error-bar range. From Fig. [2] (a-b), 
we find that the RPV weak phases of A'/3^A']'2i and ^as^^iu are not much constrained, but the 
modulus of the two RPV coupling products can be tightly upper limited. Since the SM pre- 
diction ranges of Af^p (B) in Bs K^n^ decay summarized in Table II is a little smaller 
(larger) than the corresponding measurements within la by CDF shown in Eq. ([T]), the al- 
lowed ranges of A'/32A'/]^2 and y*i2,X'iii are strongly restricted by these experimental data. We 
obtain |A;'^2^ii2I ^ [0.22, 4.86] x lO^^ and its phase (pRPv G [80°, 123°]. However, we could 
not find the allowed space of A'^^gA^^^ within la error-bar of the experimental bounds. Within 
2a error-bar of the experimental data, one can find the allowed spaces of these two RPV cou- 
pling products which are given in Fig. [2] (c-d). One can find that the RPV weak phases only 
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Figure 2: The allowed parameter spaces for the relevant RPV couplings constrained by Bs 
K~K^ and K~ii^ . (pRPv denotes the RPV weak phase. 

Table IV: Bounds on the relevant RPV couplings by Bg — > K~K~^, K~tt^ decays for 100 

GeV sfermions. Previous bounds are listed for comparison. 



Couplings 


Our bounds [Process] 


Bounds I [Process] [29] 


Bounds II [Process] [5^ 


l-'*131-*'l2ll 
l-'*'il3-*'il2l 


< 7.01 X 10--' [Bs K+K-] 

< 2.84 X 10--* [Bs K+K-] 


[0.61, 4.61x10-^ rr, j^i 
[5.6, 7.2] X 10-3 IB-*^K] 

[0.36, 1.1] X 10-3 [B ttK] 


< 1.54 X 10-2 [Bu K-TT°] 

< 2.71 X 10-3 [Bu A'-tt"] 


I^132'^112l 


< 5.01 X 10-3 [Bs K+TT-] 
[0.67, 1.90] X 10-3 [Bs K+TT-] 


[0.54, 2.9] X 10-3 [Ba tttt] 
[0.27, 0.77] X 10-3 [Ba tttt] 


< 4.69 X 10-3 [Ba ~* Tr+vr"] 

< 1.90 X 10-3 [Ba TT-TT + ] 



have the positive values, the RPV weak phase of A'/32A'/]^2 li^s in [60°, 139°] and the phase of 
X'*i-^X'iii lies in [104°, 158°]. Furthermore, the strengths of the two RPV coupling products are 
restricted strongly, which are summarized in Table IV. For comparison, the existing bounds on 



;hese quadric coupling products, which obtain from Bu,d decays o 



29 



52j ] are also listed. Note that, previous bounds-I of Ref. 29| are obtained by considering 



the same quark level process 



the experimental constraints of all relevant decay modes at the same time, so the allowed RPV 
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coupling spaces are very narrow. In Ref. 5^, the bounds are given through experimental 



restraints mode by mode. Our bounds on A'/^j^A'/gi, ^'ii?,Ki2 Xil2'^'ii2 consistent with the 



existing ones in Refs. 521], and just a little weaker than these in Ref. 29|] which are obtained 



from many correlated experimental constraints. Our bound of lA^'^gA'jj^^l also consists with one 



from Ref. 52|], however, there is only very narrow overlap between range [0.27,0.77] x 10~^ 
in Ref. |29| and ours [0.67,1.90] x 10"^, therefore, it should be of order 10~^ if jA-^gA-^J can 
survive. 

Next, we will use the constrained parameter spaces from Bg K~K^ and K~7r~^ decays, 
as shown in Fig. [2], to predict the RPV effects on the other quantities which have not been 
measured yet in Bg K^*^~ K^*^^ , K^*^~tx'^ and K^*^~ decays. With the expressions for B, 
A'^p, A^p and fi, we perform a scan through the input parameters and the new constrained 
RPV coupling spaces, and then the allowed ranges for B, A^p, A^p and /l are obtained with 
different RPV couplings, which satisfy relevant experimental constraints of Bs decays given 
in Eq. ([I]). The numerical results for B,, K^*)- K^*)+ and B^ K^*^' 71+ , K^*^' p+ are 
summarized in Table V and Table VI, respectively. 

Comparing the RPV SUSY predictions given in Table V and Table VI to the SM values 
listed in Table II and Table III, we give some remarks on the numerical results. 

• All branching ratios can be greatly changed by the RPV couplings compared to the SM 
expectations. 

• The RPV effects on A'^SpiBs ^ K*'ti+) and A'^Sp^Bs K'p^) are found to be very 
small, but could be large for the direct CPA in other five Bg K*^ p~^ and J'i"(*)~i^'(*)+ 
decays. 

• The mixin g-induced CPA in K^*)-K^*)+ system can be greatly enhanced by the 
RPV couphngs A;'*,A;'2i and X^,)^,,. 

• The squark exchange couplings A'/g^A'/gi and A'/ggA'/j^g could have significant impacts on 
fhi^Bs K*~ K*+) and fi^Bs K*~ p'^), which could be decreased as low as 0.30 and 
0.42, respectively. 

In Figs. [SB, we present correlations between the physical observable B, A'^p, A^p, /l 
and the parameter spaces of different RPV couplings by these three-dimensional scatter plots. 
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Table V: The theoretical predictions of 5, R^*)- K^*)+ for B (in units of 10"^), A\ 
and fi with the allowed regions of the different RPV couplings. 





^131^^121 


'^jl2'^jl3 


B{Bs K*~K+) 
B{Bs K~K*+) 
B{Bs K*-K*+) 


[0.33, 10.56] 
[0.01, 28.88] 
[0.29, 27.23] 


[0.32, 12.89] 


A'^'piBs K*~K+) 
A'^(i'piBs ^ K~K*+) 
A^'p'^'^iB, ^ K*~K*+) 


[-0.50, 0.50] 
[-0.25, 0.50] 
[-0.98, 0.95] 
[-0.30, 0.52] 


[-0.25, 0.23] 
[-0.11, 0.22] 


A'SpiBs ^K^K+) 
A'gfiBs^Bs K*-K+) 
A'SpiBs&iBs K'K*+) 
A'S^'''{Bs ^ K*-K*+) 


[-0.97, 0.98] 
[-0.89, 0.97] 
[-0.89, 0.96] 
[-1.00, 1.00] 


[-0.99, 1.00] 
[-0.98, 0.55] 
[-0.98, 0.58] 


fUBs ^ K*-K*+) 


[0.30, 0.97] 





Table VI: The theoretical predictions of K^*^'7r+ , K^*^' p+ for B (in units of 10"^), 

Aqp and fi with the allowed regions of the different RPV couplings. 





^132^112 


^il3'^ill 


B{Bs K*~7r+) 


[1.34, 11.59] 


[4.49, 13.47] 


B{Bs ^ K~p+) 


[2.02, 23.32] 




B{Bs ^ K*~p+) 


[0.40, 3.31] 




Af^UBs K*~7r+) 


[-0.04, 0.04] 


[0.00, 0.02] 


Ai^^Bs - K~p+) 


[-0.05, 0.01] 




A^c'/'"iBs - K*-p+) 


[-0.25, 0.73] 




fUBs - K*~p+) 


[0.42, 0.96] 
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From Figs. [30, one can see the changing trends of the physical observables with the modulus 
and RPV weak phase (pRpv- Taking the first plot in Fig. [3]^a) as an example, this plot 
shows B{Bs K*~ K^) changing trend with RPV coupling A'/g^A'/gi, where projections on 
three perpendicular planes are also given. The jA'/g^A'/gil-^/jpy plane displays the allowed 
regions of A']'3]^A']'2i which satisfy experimental data in Eq. ([T]) (the same as the Figl2](a)). The 
B{Bs /i*-i^'+)-|A;'^iA'/2il plane shows that B{Bs K*- K+) is increasing with IX'isiXuil, 




Figure 3: The effects of RPV coupling A'/^^A'/ai in R- K+ , K' K*+ , K*' K+ , K*-K*+ decays. 

B in units of 10~^ and lA'/oi A'/ni I in units of 10^^. 
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"4/ " 
(d) 

Figure 4: The effects of RPV coupling Af^gA'^is in Bs K~ K+ , K*+ decays. B in units of lO'^ 
and lAj^gA'j^gl ™ units of 10"'^. 

the B{Bs K*-K+)-^ 

RPV plane shows that B{Bs ^ K* K^) is decreasing with |0_Rpy|. 
Additional refined measurements of B{Bs K~ K^) can further restrict the constrained space 
of A']'3]^A'/2i, thus more accurate B{Bs K*~ K^) can be predicted. The following salient 
features in Figs. OEl are summarized as following. 

• Fig l3] displays the effects of RPV coupling A'/g^^A'/gi on B, A^p, A^p in penguin-dominated 

Bs K-K+,K-K*+, K*-K+ and B, A^p\ A^p'"' , /l in penguin-dominated 5, ^ 
K*^ K*^ decays. The constrained |A'/3]^A'/2i|-0Rpy plane shows the allowed range of 
^131^121 as shown in Fig. [2](a). The ^ K*~K+, K-K*+, K*-K*+) shown in 

Fig. [3](a-c) have the similar change trends with |A'/3^A'/2i| and l^ppy], and they all in- 
creases with |A'/3;^A'/2il and decreases with \(pRpv\- For the Aqp/Aqp^ shown in Fig. [3](d- 
g), A'/3^A']'2i coupling contribution could be significant. lAf^pi^Bg K~K~^)\ increases 
when |A'/3]^A']'2il is small, then Aqp{Bs — ^ K~ K^) decreases and its sign is flipped. 
A'}^'p{Bs(Bs) K*-K+), Afp{Bs(Bs) ^ K~K*+) and Acp\Bs K*~ K*+) could 
have smaller range with |A'/^iA;'2i|. |^cp(^^ ^ K"K+,K*"K+)\ and \A^p\Bs ^ 
K*-K*+)\ decrease with |0ppy|. A'^Sp^Bs ^ K*-K+,K-K*+) and^^p''(^« K*~ K*+) 
could be close to zero in entire (j)Rpv range. As shown in Fig. [3](h-k), four mixing-induced 
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(g) 



Figure 5: The effects of RPV coupling A']'32A']'i2 K* 7r+, K^p+, K*~ decays. B in units 

of 10"^ and lA'/SoA'/iol in units of IQ-^. 




(a) (b) 



Figure 6: The effects of RPV couphng A ■'J^^ A ■5^3 in B ^ K* 7r+ decays. B in units of 10 ^ and 
l-^iis-^iiil units of 10"^. 



CP asymmetries are very sensitive to |0_Rpy| but not sensitive to |A'/3]^A'/2il- For the pen- 
guin dominated process Bg K*~K*~^, its longitudinal polarization could be small as 
shown in Fig. M}), however, most points of K*~ K*^) fill in [0.7,0.9]. 
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• The effects of Af^gA^^a on B, A'^Sp and ^^)? of ^ R- K+ , K*- K+ , R- K*+ are pre- 
sented in Fig. m Tlie constrained \^'i2?,Ki2\''t'RPV plane is tlie same as Fig. [2]^b). Fig. 
m^a) sliow tliat B{Bs K*~ K'^) increases witli lA'j^gA^^^al ^^d decreases witli \(t)Rpv\- 
As sliown in Fig. Hl^b), at first \Acp{Bs K" K^)\ increase with lA-^JgA-^gl' ^^en 
Acp{,Bs K~K^) could occupy the entire range [-0.25,0.23] when |A-i3A-i2l lies in 
[1.6,2.8] X 10-3, and \A'^Sp{Bs ^ K- K+)\ decreases with |0/jpy|- Afp{Bs K*- K+) 
has narrow ranges with the constrained |A'j'^3A^]^2l ^nd jf/'ijpyl- Fig. Hl^d-f) show the RPV 
effects in relevant mixing-induced CPA. \Aqp{Bs K~K~^)\ is not sensitive to |A'-'^3A^-|^2l 
but sensitive to ^rpv- lA'SfiB.LBs K*~ K+)\ and ^ K-K*+)\ de- 
crease with lA'j'^gA'jial ^^d they first increase and then decrease with |0i?py|. 

• In Fig. El we plot B, Af}'p of B,, A'*-7r+, K~p+, and B, Acf", fi of 5, ^ K*- p+ 
decays as functions of A']'32Aii2- The constrained |A'/32A"]^2l"</'-RPV plane is the same as 
Fig. [2]^c). One can find B{Bs /i*"7r+, increase with |A'/32A'/i2l ^^d \(t)Rpv\- 
B{Bs K*^ p^) first decreases and then increases with |A'/32A'/;^2l5 ^^d it is not very 
sensitive to |0/jpy|. As shown by Fig. El^d-e), the squark exchange RPV effects on 
Acp{B K*~'7i~^ , K~ p'^) are very small. Acp{B K*~'7i~^ , K~ p'^) first decrease and 
then increase with |A'/32A']'^2l) ^^d they both increase with (pppv- A^'p^{B — > K*~ p^) is 
sensitive to A"32A']'i2 coupling, and could be enhanced to ~ 70% when |A']'32A']'i2l is around 
3 X 10"'^. Aqp^{B K*~ p^) first increases and then decreases with |A'/32A']';^2l5 but is 
not sensitive to (pppv- fii^B K*~p~^) has the largest allowed range when |A'/|2Aii2l is 
around 3 x 10"^. The A'/32A'/i2 couplings could decrease — > K*~p~^) to 0.42. 

• Fig. E] shows the effects of the RPV couplings X'*^^^^-^^^ in Bg K*~7r^ decay. B{Bs 
K*~'n'~^) increases with |A^*3A^]^;^| and is insensitive to (pppv- A'^jpiB K*~Tx+) decreases 
with |A^*3A'ji]^| and increases with (pppv- 

4 Conclusions 

In conclusion, we have studied the eight decay modes Bs (*)+^ K^*'>~7r~^, K^*^~ p'^ 

in the RPV SUSY with the QCDF for the hadronic dynamics. With the recent experimental 
data of Bg decays, we have obtained fairly constrained parameter spaces of the RPV couplings. 
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Furthermore, using the constrained parameter spaces, we have shown the RPV SUSY expec- 
tations for the other quantities in Bg K^*^~ K^*^~^ , K^*^~7r~^, K^*^~ decays which have not 
been measured yet. 

We have found that the RPV couphngs A'/3^A']'2i and \'*i-^\'ii2 could significantly affect 
penguin-dominated Bs decays. Within the parameter spaces already highly 

constrained by Bg K~K~^, the branching ratios of B^ K^K*^ , K*~K^ and K*~ K*~^ 
could be enhanced by few times, and the direct CPA and the mixing- induced CPA are in quite 
large ranges. Interestingly, the longitudinal polarization fraction of Bg K*~ K*^ could be 
suppressed as low as 0.30. Therefore future experimental measurements of these decays could 
shrink or reveal the relevant NP parameter spaces. It is found that the squark exchange cou- 
pling A'^'gjA'/^g could have large contributions to the branching ratios of Bg /i*"7r+, p~^, 
and enhance the longitudinal direct CP asymmetry of Bs K*~ to ~ 70%. The longi- 
tudinal polarization fraction of Bs K*^ p^ could be suppressed too. The slepton exchange 
coupling X'*i^Xiii could enhance the branching ratio of Bg K*~tt^ by few times. We also 
have presented correlations between these physical observable quantities and the constrained 
parameter spaces of RPV couplings in Figs. [3H6l The results in this paper could be useful for 
probing RPV SUSY effects and searching direct RPV signals at Tevatron and LHC in the near 
future. 
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